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'""' ' Abstract 

X5 ' 

D . We study the dynamics of anisotropic Bianchi type-IX models with matter and 

1^ I cosmological constant. The models can be thought as describing the role of anisotropy 

in the early stages of inflation, where the cosmological constant A plays the role of 

the vacuum energy of the inflaton field. The concurrence of the cosmological constant 

j—{ ', and anisotropy are sufficient to produce a chaotic dynamics in the gravitational de- 

K^ I grees of freedom, connected to the presence of a critical point of saddle-center type in 

^T ' the phase space of the system. In the neighborhood of the saddle-center, the phase 

(^ , space presents the structure of cylinders emanating from unstable periodic orbits. The 

Cn I non-integrability of the system implies that the extension of the cylinders away from 

1^ ■ this neighborhood has a complicated structure arising from their transversal cross- 

Q^ . ings, resulting in a chaotic dynamics. The invariant character of chaos is guaranteed 

Q I by the topology of cylinders. The model also presents a strong asymptotic de Sitter 

O^' attractor but the way out from the initial singularity to the inflationary phase is com- 

^ . pletely chaotic. For a large set of initial conditions, even with very small anisotropy, 

. . I the gravitational degrees of freedom oscillate a long time in the neighborhood of the 

, ^ ' saddle-center before recollapsing or escaping to the de Sitter phase. These oscillations 

^ ■ may provide a resonance mechanism for amplification of specific wavelengths of inho- 

H , mogeneous fluctuations in the models. A geometrical interpretation is given for Wald's 

■ " " ' inequality in terms of invariant tori and their destruction by increasing values of the 

cosmological constant. 



1 Introduction 

One of the cornerstones in the paradigm of inflation [|l| is the presence of the cosmo logical 
constant, arising as the vacuum energy of the inflaton field. The cosmological constant plays 
a fundamental role in the gravitational dynamics of the inflationary model, by inducing an 
exponential expansion of the scales of the model toward the de Sitter configuration. This 
asymptotic approach to the de Sitter solution constitutes the basis of the so-called cosmic no- 
hair conjecture. In the realm of homogeneous cosmologies, Wald0 showed that all initially 
expanding Bianchi cosmological models with positive cosmological constant, except Bianchi 
IX, evolve towards the de Sitter configuration. Bianchi IX models demand further that 
the absolute value of the cosmological constant be sufficiently large compared with spatial 
curvature terms. For the more general case of inhomogeneous models, Starobinskiip] showed 
that they do inflate if a positive cosmological constant is present. This crucial aspect of the 
cosmological constant has been sufficiently emphasized in the literature, and many authors 
have examined its role in producing non-trivial dynamics in the early stages of inflation. 
In particular, Calzetta and El Hasi[^, and Cornish and Levin|p exhibited chaotic behavior 
in the dynamics of Friedmann-Robertson- Walker models with a cosmological constant term 
and scalar fields conformally and/or minimally coupled to the curvature. Due to this feature 
of the dynamics, small fluctuations in initial conditions of the model preclude or induce the 
universe to inflate. The Hamiltonian dynamical system originating from the fleld equations 
is complex. In Ref. [4], the effective degrees of freedom of the model are the scale factor and 
a conformally coupled scalar fleld. Such a scalar fleld is interpreted as a radiation fleld, which 
is later assumed to gain mass by the presence of the inflaton fleld. In Cornish's work, the 
degrees of freedom are the scale factor and two minimally and/or conformally coupled scalar 
flelds. In both cases, the cosmological constant is responsible for the existence of a critical 
point S in the flnite region of the phase space of the model. S" is a pure saddle point [|| and the 
separatrices emanating from S are wholly contained in an invariant plane of the dynamical 
system, corresponding to the gravitational degree of freedom only. The separatrices connect 
S to other critical points, and are denoted heteroclinic|TD|. These connections are known 
to be highly unstable, and their breaking, due to the perturbations come from the coupling 
of the gravitational variable with the scalar flelds, are basically responsible for the chaotic 
dynamics referred to above. This chaotic behavior is consequence of the so-called Poincare 
homoclinic phenomena in dynamical systems [p]]. 

There is however another important feature in the pre-inflationary dynamics, arising 
from the presence of a positive cosmological constant, whenever anisotropy is also present 
even in the form of small perturbations. This point has not been emphasized yet in the 
literature of inflation and will be object of our paper. The degrees of freedom of the system 
are taken basically in the gravitational sector and, for simplicity we restrict ourselves to 
two distinct scale factors. The conjunction of the cosmological constant and anisotropy 
implies the existence of a critical point E in phase space, identifled as a saddle-center |^. 
As a consequence, we have a wealthy dynamics based on structures as homoclinic cylinders 
which emanate from unstable periodic orbits that exist in a neighborhood of E. Analogous 
to the breaking and crossing of heteroclinic curves, for instance, these cylinders will cross 
each other in non-integrable cases producing a chaotic dynamics. As we shall discuss, these 
structures will constitute an invariant characterization of chaos |T^ in the models. We will 



also discuss their implications for the occurrence or not of inflation, as well as for the physics 



in the early stages of inflation. The paper is organized as follows. In Sec. II we establish 
the Hamiltonian and the basic characteristics of the model. Sec. Ill is devoted to present 
the cylindrical topology near the saddle-center which will be very useful to understand the 
dynamics of orbits on the phase space. The numerical evidence for a physically relevant 
chaotic behavior is showed in Sec. IV, whereas in Sec. V an interesting connection between 
the break up of torus and the Wald's analysis is discussed. Finally, in Sec. VI, we conclude 
and trace some perspectives of the present work. 

2 The Dynamics of the Model 

We consider anisotropic Bianchi IX cosmological models characterized by two scale functions 
A{t) and B{t) with the line element 

ds^ = dt' - A\t) {w'y - B\t) [{w^ + {w^)% (1) 

Here t is the cosmological time and {w^, w"^, w^) are invariant 1-forms for Bianchi IX models|T 



The matter content is a assumed to be a perfect fluid with velocity field 6q in the comoving 
coordinate system used plus a cosmological constant A. The cosmological constant term 
is interpreted as arising from the vacuum energy of the inflaton field such that our mod- 
els may provide a description for a pre-infiationary anisotropic stage of the universe. The 
energy-momentum tensor of the fluid is described by 

T>^'^ = {p + p)6i;6^o-P9'"' (2) 

where p and p are the energy density and pressure, respectively. We assume the equation 
of state p = 7P, 0<7<1. For sake of simplicity, we restrict ourselves to the case of dust 
(7 = 0). Distinct features arising in the cases of 7 7^ will be discussed in the conclusions. 
Einstein's field equations [p^j F] 



G^^-A^'^^ = T^^ (3) 

for Eqs. (1) and (2) can be obtained from the Hamiltonian constraint 
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HiA, B, Pa, Pb) = -^^ - —^ + 2 A - -— - 2 AAB' - Eo = 0. (4) 



where Pa and Pb are the momenta canonically conjugated to A and B, respectively, and Eq 
is a constant proportional to the total energy of the models. It also occurs as the first integral 
of the Bianchi identities, 2pAB'^ = Eq. The full dynamics is governed by the Hamilton's 
equations 
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The dynamical system (5) has one critical point E in the finite region of the phase space 
whose coordinates are 

E: Ao = Bq = ^;Pa = Pb = 0, (6) 

V4A 

with associated energy Eq = Ecr = ~^- This critical point represents the static Einstein 
universe. Linearizing (5) about the critical point E, we can show that the constant matrix 
determining the linear system about E has the four eigenvalues 

Ai,2 = ±^;A3,4 = ±^, (7) 

which characterizes the E as a saddle-center. The system (5) has also a degenerate critical 
point [l^ at A = B = 0, Pa = Pb = 0. A straightforward analysis at the infinity of the phase 



space under consideration shows that it has two critical points at this region, corresponding 
to the de Sitter solution, one acting as an attractor (stable de Sitter configuration) and the 
other as a repeller (unstable de Sitter configuration). The scale factors A and B approach to 
the de Sitter attractor as A = S ~ eV^* and Ps = 2 Pa ^ e2^/A/3t_ q^^^ ^^ ^^^e questions 
to be examined in this paper is the characterization of sets of initial conditions for which 
this asymptotic de Sitter attractor is attained. 

Another important feature of the dynamical system (5) is to admit an invariant manifold 
Ai defined by 

A = B, Pa = Pb/2. (8) 

On Ai the dynamics is governed by the two-dimensional system 

Pb 



B = B^ 



8^ ^ (9) 



I Pb = 4w-1 + ^AB^ 

The system (9) is integrable. Its phase portrait is depicted Fig. 1, where the integral curves 
represent homogeneous and isotropic universes with Bianchi IX compact spatial sections, 
with dust and a cosmological constant. 

We remark that the critical point E belongs to Ai. Contrary to the models examined in 
Refs. Q and P , the separatrices S are present as solutions of the full dynamical equations. 
As we will discuss in Sec. 3 and 4, the behavior of the separatrices found in the previous 
studies will be played by the cylinders emanating from unstable periodic orbits associated 
to the saddle-center of E. 

3 Cylindrical Topology near the Critical Point E 

As a consequence of the saddle-center nature of the critical point E, the dynamics of the 
models described by (5) exhibits completely new features, which will be analyzed now. Part 
of this analysis is based on Refs. 0||3][§1- The overall scenario originates from both the 
anisotropy of the model and the cosmological constant. To discuss the topology of the phase 
space in the neighborhood of E, we make use of a theorem by Moser|T^ which states that 



it is always possible to find a set of canonical variables such that, in a small neighborhood 
of a saddle-center at the origin, the Hamiltonian is expressed as 



H{q,,q2,PuP2) = ^{pI- ql) - v^(p? + ql) + C(3) - ^o + E,, = 0. (10) 

The critical point is located at the origin qi = q2 = Pi = Pi = 0, with Eq = Ecr- Here C(3) 
denotes non-quadratic terms of the expansion and ± vA, ± 2 z v2A are the eigenvalues of 
the linearized system about E. 

Let us now restrict ourselves to a small neighborhood of E such that we can neglected 
0{3) and the quantity S = Eq — E^r is small. The Hamiltonian reduces to 



Hiqi, q2,Pi,P2) ~ ^ (P2 - q!) -V2Aipl + ql) - ^ = 0. (11) 

In this approximation, H is separable with two approximate constants of motion given by 
the partial energies 

E2 = ^ {pI - ql), E, = v/=2A (pI + ql) (12) 



The energies E2 and Ei will be referred ||T7[| to as the hyperbolic motion energy and the 
rotational motion energy of the system about E, respectively. Note that Ei is always positive. 
To describe all possible motions, the following situations must be taken into account. 

If £^2 = two possibilities arise. First, we have P2 = Q'2 = meaning that the motions 
are unstable periodic orbits teq in the plane (pi, qi). Such orbits depend continuously on the 
parameter Ei ~ —£ (cf. Fig. 2(a)). The second possibility will be p2 = ±(3'2, which defines 
the linear stable Vg and unstable Vu 1-dimensional manifolds of Fig. 2(b). These manifolds 
are tangent at E to the separatrices S of the invariant manifold (described as gi = 0, pi = 
in the new variables) of Fig. 1. The separatrices are actually the non- linear extension of 
Vu and Vg. The direct product of the periodic orbit te^ with K, and Vu generates, in the 
linear neighborhood of £^, the structure of stable {teq x Vs) and unstable cylinders {teq x Vu), 
which coalesce into the orbit teq for times going to +cxd or — cxd, respectively (cf. Figs. 2(c), 
2(d)). The energy of any orbit on these cylinders is the same of the periodic orbit teq- It 
can be showed that, in the non-linear regime (when non-quadratic terms of the Hamiltonian 
must be taken into account), the plane-(gi,pi) of the rotational motion in the linear regime 
extends to a 2-dimensional manifold - the center manifold^ - of unstable periodic orbits of 
the system. The intersection of the center manifold with the energy surface 

H{S) = (13) 

(cf Eq. (11)) is a periodic orbit parametrized with Eq, from which a pair of cylinders emanate. 
We can see easily from Eq. (11) that the intersection of the central manifold with the energy 
surface £^ = is just the point E. For S > the energy surface does not intersect the central 
manifold, the occurrence of the structure of cylinders being therefore restricted to the energy 
surfaces in which S < 0. 

For the case E2 y^ and S < 0, the motion is restricted on infinite cylinders resulting 
from the direct product of the hyperbolae lying in the regions / and // of Fig. 2(b), with 
periodic orbits of the central manifold in a small neighborhood of E. A general orbit which 



visits the neighborhood of E belongs to the general case Ei ^ E2 ^ 0. In this region the 
orbit has an oscillatory approach to the linear cylinders (cf. Fig. 3), the closer as E2 — > 0. 
For instance, the outcome of this oscillatory regime will be collapse if i^^2 < or escape to 
de Sitter attractor if ii^2 > for initial conditions taken in the quadrant ^2 < 0, P2 > 0. 
This is for the linear regime. In general, for orbits which visit a neighbohood of E, the non- 
integrability of the Hamiltonian system (5) induces that the partition of the energy S into 
the rotational mode energy Ei and the hyperbolic mode energy E2 is chaotic. In another 
words, given a general initial condition of energy S, we are no longer able to foretell which 
amount of S goes to each mode, namely, in what of the regions / or // (Fig. 2(b)) the 
orbit will land when approaches of E. This manifestation of the non-integrability will be 
physically meaningful to characterizing a chaotic exit to inflation, as we will show in the 
next section. 

We must finally discuss the extension of the structure of cylinders outside the neighbor- 
hood of E. Let us consider a linearized cyhnder associated to a periodic orbit with energy 
Eq such that £ = Eq — Ecr is very small. One possible way to examine the extension of this 
cylinder is to consider the linearization of the dynamical system (5) about the separatrix S. 
The equations for the separatrix are denoted by 



{as{t),ps{t)), 



(14) 



and consider the expansion about S 



as + X, 
= Ps + Z 

as + Y, 
= 2ps + W. 

To the first order in (X, Y, Z, W), the system (5) results 




(15) 
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where A(t) is a 4 x 4 time-dependent matrix, the entries of which are simple functions of 
as and ps- For the sake of completeness we give its expression in the Appendix, but for 
our purposes here it will be enough to remark that A(t) is bounded, except for 05 — *> 0. A 
matrix is defined to be bounded is all its entries as well as its determinant are bounded^. 

Clearly, for a small neighborhood of E, the separatrix is approximated by Vs and Vu, and 
the linear cylinder based on teq will be a solution of (16) by construction |jl8|. Furthermore, it 
follows from the system (16) that its extension will be contained in a 4-dimensional small tube 
about the separatrix, as far as as(t) does not go to zero. When as(t) — > the above linear 
approximation is no longer valid. Higher order terms become important for the dynamics, 
and the non-integrability of the system results in the distortion and twisting of the cylinders. 



^A similar analysis can be made of the motion around other integral curves of the invariant manifold 
depicted in Fig. 1. In Eqs. (14) and (15) it is sufficient to substitute as{t) and ps{t) by the solution of the 
integral curve (aj{t),pj{t)) associated to a given energy Ej in the invariant manifold. 



The stable cylinder and the unstable one will cross each other, producing chaotic sets00 
and consequent "destroyed" regions of the Poincare maps of the system, as showed in Section 
5. This also occurs for cylinders emanating from periodic orbits of the center manifold which 
are not in a small neighborhood of E (periodic orbits in the non-linear regime). 

Finally, with view to the next section, let us select a 4-dimensional small sphere of initial 
conditions about one point 5*0 of the separatrix with radius R{So) of the order of the linear 
perturbation in (15), namely, R{So) ~ (X^ + Y^ + Z^ + W'^fl'^ (cf. Fig. 4). The energy 
surfaces which intersects the sphere are those with Eq in the domain |ii^o~-^cr| < R{So)- The 
cylinders (associated to the periodic orbits te^) in these surfaces will obviously intersect the 
sphere. As we discussed before, if we evolve the sphere back as ^ the small tube spread 
and twist, as well as the cylinder associated to teq- This stable cylinder and the unstable 
one will cross each other and return eventually several times to the sphere. This geometry 
will be the basis of the numerical experiments of the next section. By evolving dynamically 
the sphere towards the neighborhood of E, we show that for a band of energy, the sphere 
contains a chaotic set which induces a chaotic exit to inflation. 



4 Chaotic Exit to Inflation 

In the numerical experiments performed here, we use the variables the variables {A, B, Pa, Pb) 



All calculations were made using the package Pomcare|jT9| in a IBM compatible PC Pentium 
133 with 64 MB of RAM memory, and a Fortran program to construct Poincare surface of 
sections, where we enforce that the error of the Hamiltonian never exceed a given threshold 
of 10"^". Henceforth, we assume A = 0.25, so that the critical point E is characterized by 
A = B = 1.0, Pa = Pb = Q and E„. = 1. 

he phase space under consideration is not compact, and we will actually identify a chaotic 
behavior associated to the possible asymptotic outcomes of the orbits in this phase space, 
namely, escape to de Sitter state attractor at infinity (inflationary regime) or collapse after a 
burst of initial expansion. To begin with our numerical experiments, we consider the invariant 
manifold 7V1 (A = i?, P4 = Pb/2). The basic characteristic of the curves in M. have been 
discussed already (cf. Fig. 1). We remark, however, that the separatrices define the regions 
of collapse and expansion in M., but not in the full 4-dimensional phase space. Following 
the theoretical backgrounds presented in Sec. 3, we investigate numerically the behavior 
of orbits in a domain near the separatrices. In essence, the sets of initial conditions are 
constructed in the following way. Let 5*0 be a point belonging to the separatrix [Eq = 1.0), 
with coordinates A = B = 0.4, P^ = 2 Pa = 1.357645019878171. We are interested, 
therefore, in those orbits representing expanding models after the initial singularity. Around 
this point, we construct a 4-dimensional sphere in the phase space with arbitrary small radius 
R, for instance R = 10~^ 10"^, 10"^, etc (cf. Fig. 4). The values of A, B, Pa and Pb are 
taken in energy surfaces which have a non-empty intersection with this sphere as evaluated 
from the Hamiltonian constraint. As we have seen, such energy surfaces are those for which 
the range of energy, A Eq, about Eq = Ecr = 1.0 is of the order of, or smaller than the radius 
R. On physical grounds, we are considering the evolution of cosmological models with small 
anisotropic perturbations. 

After several numerical experiments with the above sets of initial conditions, we note that, 
as expected, two possible outcomes arise: collapse or expansion into de Sitter configuration. 



The energy Eq that varies from 1.0 — A £"0 to 1.0 + Ai^^o, with AEq ~ i?, determines the 
long-time behavior of the orbits. In Figs. (5) and (6), collapse and escape of 50 orbits are 
displayed for R = 10~^. The orbits oscillate around the separatrix and approach the critical 
point E inside a sphere about this point with radius of the same order of R. In this region, 
the linear approximation is valid and the local separation of the dynamics into rotational 
motion and hyperbolic motion (cf. Eq. 10) can be used to understand the results. The 
final state of the orbits depends crucially on the partition of the total energy S = Eq — E^r 
into the rotational motion mode and the hyperbolic motion mode. Hence, if ii^2 > the 
orbits escape, whereas collapse is characterized by E2 < 0. The rotational motion (energy 
El) describes the oscillatory character of the orbits around the critical point, indicating that 
some orbits can spend more time around E than others. 

The main result of this section is to show that, for a determined interval of energy 6 E* 
contained in the domain V = [1.0 — AEq, 1.0 + AEq], this partition is chaotic. Indeed, let 
us try to determine in V the values of energy Eq = Emm for which all orbits escape and 
£^0 = Emax for which all orbits collapse. According to our numerical work, we find out that for 
each sphere of initial conditions, there always exists a non-null interval S E* = \Emax — Emin\, 
where part of the orbits escapes and another part collapses, resulting in an indeterminate 
outcome. In Fig. (7), this behavior is showed for spheres of initial conditions of radius 
R = 10~^, 10~^ and 10~^. An empirical relation between the gap 6 E* and the radius R 
is obtained and showed in Table 1. For spheres with R < 10^^, we have S E* (x R^. 
Indeterminate outcome due to 6 E* y^ occurs only for S = Eq — Ecr < 0, as expected 
(indeed, if £^ > all orbits escape since, from Eq. (11), E2 > 0). The orbits collapse or 
escape, depending on the partition of the energy £ into the modes Ei and E2 (cf. Eqs. 11 
and 12), such that the later assumes negatives or posititves values, respectively. The way 
in which this partition works, once a set of orbits approach E, is completely unknown for 
energies in within the gap 6 E*. As a consequence, any infinitesimal fluctuation from a given 
initial condition inside the sphere can lead to an indeterminate outcome, that is, collapse or 
escape. This is the evidence of chaos in a physically relevant context. In other words, we 
may state that the boundaries of initial conditions for collapse and inflation are mixed as a 
consequence of the crossing of cylinders. 

The presence of chaos in the system is a consequence of the crossing of stable and un- 
stable cylinders, emanating from the unstable periodic orbits of the center manifold, as 
dicussed before. This topological structure is actually an invariant characterization of chaos. 
Finally, we remark that the above behavior is not restricted to initial conditions taken in 
small neighborhoods of points of the separatrix. Any sets of initial conditions taken in an 
arbitrary neighborhood of the invariant manifold Ai which result in orbits that visit a small 
neighborhood of E, display the above chaotic behavior. 

In a less simple model, radiation could be taken into account as a more plausible matter 
field emerging after the Planck era||2D| instead of dust. The effect of radiation is considered 
if we set the equation of state as p = | p. The saddle-center critical point and the invariant 
manifold are also present, and the same type of behavior occurs if we construct sets of 
initial conditions as before. There is also a similar relation between the gap 6 E* and the 
radius R showed in the Table 2. A more complete and detailed analysis will be subject of a 
forthcoming paper. 



5 Further Numerical Results and Wald's Phenomenon 

We consider here the evolution of completely anisotropic models taking initial conditions far 
from the invariant manifold Ai. It is a remarkable fact that the Hamiltonian (4) is regular 
at the plane A = where the curvature is singular. Therefore, from the point of view of 
the Hamiltonian dynamics, orbits can be analitycally extended beyond the singularity to the 
domain A < 0. Although these orbits are not physically meaningful in this domain, they 
are nevertheless essential in the description of the underlying geometrical structure of the 
full Hamiltonian system and its non-integrability. In this context we will give a geometrical 
picture for Wald's inequality0. 

In Fig. 8(a), we exhibit the Poincare map of the system for {Eq = 0.92504831163435, A = 
0.25) with surface of section (Pg = 0, Pg > 0). This map makes explicit the coexistence of 
KAM tori and destroyed regions in the phase space, as a consequence of the non-integrability 
of the system. We were not able to find a torus totally contained in the domain A > 0. Thus, 
orbits on, or inside the tori (with initial conditions taken on A > 0) necessarily collapse. 
This can be seen from double Poincare map with surface of section Pb = as showed in Fig. 
8(b), where {Eq = 0.15, A = 2.0). On the contrary, orbits in the destroyed region are free to 
escape. Orbits on a torus with analytical continuation to a domain of negative A are said to 
perform cosmic cycles (according with Ref. [Q). 

Our next step is to examine the effect of varying the value of the cosmological constant 
A on the tori structure in the same region of phase space. As can be verified numerically, 
the effective result of increasing A is the destruction of the tori. We showed this in Fig. 9, 
where orbits are drawn in the same region of initial conditions but with distinct values of 
A. If we change A from 0.25 to 0.30 (the energy changes from Eq = 0.8, taken initially, 
to Eq = 0.8000050001250895), orbits initially in a torus escape to the de Sitter attractor 
after some cycles; if we further increase A, no cosmic cycles takes place before escape. The 
phenomenon of destruction of tori, by increasing the value of the cosmological constant, is a 
geometrical picture for Wald's inequality. We nevertheless remark that, for any value of the 
cosmological constant (however large, except infinity), there will always be collapsing orbits. 

6 Final Remarks and Conclusions 

In this paper we have discussed the dynamics of Bianchi-IX models, which may provide a 
description of pre-inflationary stages of the universe. The main ingredients of the models 
are a cosmological constant, arising as the vacuum of the inflaton field, and anisotropy. 
The degrees of freedom are taken in the gravitational sector, and for simplicity we restrict 
ourselves to two distinct scale factors only. The presence of anisotropy, even in the form of 
small perturbations, together with the cosmological constant produces in the phase space 
of the system a critical point identified as a saddle-center. Associated to the saddle-center 
we have a 2-dimensional manifold of unstable periodic orbits, and the structure of infinite 
unstable cylinders emanating from them. The stable and unstable cylinders coalesce into 
the corresponding periodic orbit for time going to oo and — oo, respectively. The non- 
integrability of the system results in the distortion and twisting of the cylinders, and the 
eventual intersection of one (stable) with the other (unstable) with a consequent chaotic 
behavior of the dynamics of the phase space. This behavior of the cylinders is analogous 



to the one played by the separatrices (heterochnic curves connecting saddle points) of the 
gravitational sector in the models discussed in f^ and f^. In these models the gravitational 
variable dynamics couples with the scalar fields giving rise in the heterochnic breaking and 
tangle which is the origin of chaos in these models. 

The results of Section 4, where we describe chaotic exit to infiation, extends the result of 
Cornish and Levin ||^ for the case of two gravitational degrees of freedom. The breaking of the 
boundary between initial conditions domains of collapse and escape to infiation, showed in |^ , 
is indeed meaningful for the case of one gravitational degree of freedom only (the unperturbed 
separatrix is in fact a sharp boundary between the domains of initial conditions). In our 
case, the separatrices are present in the full dynamics and define sharply regions of collapse 
and infiation in the invariant manifold, but not in the full 4-dimensional phase space of the 
gravitational dynamics. However, for each small sphere of initial conditions taken about one 
point of the separatrix, it is always possible to find a small domain (or gap) of energy such 
that the intersection of the sphere with an energy surface associated to the above domain 
is a chaotic set in the sense of discussed in Section 4. Namely, a small perturbation in 
initial conditions taken in this set would change an orbit from collapse to escape to the de 
Sitter phase. Furthermore, we conjecture that it is feasible to select initial conditions in 
this set such that the scale factors oscillate an arbitrary fixed time T about E {T = oo 
included) before collapsing or escaping to a de Sitter phase. Putting in another words, chaos 
is established by the incertainty in the partition of the energy Eq into the rotational motion 
energy and hyperbolic motion energy about E. 

The oscillations of the scale factors about E may have an important physical effect 
concerning inhomogeneous perturbations. Let them be scalar field perturbations and/or 
matter density perturbations in this gravitational background. In fact, the time- dependent 
amplitude of each Fourier component of the perturbation will satisfy a linear differential 
equation, the coefficients of which will be oscillatory functions about E. As we have seen, 
the latter may be approximated by a periodic function describing teq- Therefore, by a 
resonance mechanism, there will occur amplification of the particular Fourier components 
having period approximately equal to the period te^ (we obviously assume here that no 
dissipation effects are included in the equations of motion of the perturbation in this stage). 
Even if the universe inflates afterwards, the relative ratio of amplitudes produced after this 
mechanism of amplification will be nevertheless mantained. 

Finally, as discussed in Section 4, we have given a geometric interpretation for Wald's 
result in terms of invariant tori and their destruction by increasing of the cosmological 
constant. For a given value of the cosmological constant, the phase space presents the 
structure of KAM tori, and orbits on, or inside these tori collapse necessarilly. However, if we 
increase the value of the cosmological constant, the tori are destroyed and orbits previously 
on, or inside them may eventulally escape to de Sitter configuration. 
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8 Appendix 

The matrix A{t) is given by 



A 



4ag A as 

4 -ii + Sr + ^^^s) fir, 
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Figure Captions 

Fig. 1 Integral curves on the invariant manifold A = B, Pa = Pb/'^- The orbits on regions 
(/) have Eq < -7=^, whereas for those orbits on (//), Eq > —^- The separatrices are 
characterized by the energy Eq = E^r = /jj^ . 

Fig. 2 (a) Periodic orbit of the Hamiltonian system (5) in the linear approximation and 
projected onto plane-(gi,pi) of the normal variables, (b) The linear unstable Vu and stable 
Vs 1-dimensional manifolds. The hyperbolae are the linearized solutions in the plane (i?, Pg) 
of the saddle for E2 <Q (region /) and £^2 > (region II) . (c) Stable and unstable cylinders 
manifolds emanating from the periodic orbit r. They are the non-linear extension of the 
linearized cylinders t ^Vu and r x V^ in the neighborhood of E. (d) Numerical ilustration of 
the linear stable and unstable cylinders for A = 0.25 and Eq = 0.9999999 in the neighborhood 
oiE. 

Fig. 3 General orbits with oscillatory approach to the cylinders on the neighborhood of 
E corresponding to Eq = 0.99999999. on the plane (92,^2) of the above orbits. Note the 
asymptotes \4 and Vs. 

Fig. 4 Projection of the 4-dimensional sphere of initial conditions about the point So of 
the separatrix S. 

Fig. 5 Collapse of 50 orbits initially inside an sphere of initial conditions of radius 
R = 10"^ for Eo = 0.9999999937 and A = 0.25. (a) View of orbits projected onto the 
plane- {A, Pa), (b) Zoom of the region near the critical point. Note the oscillations around 
the separatrix as well as the critical point E. 

Fig. 6 Escape of 50 orbits to the de Sitter configuration. We consider an sphere of initial 
conditions of radius R = 10~^ within an energy surface Eq = 0.99999999999. (a) View of 
orbits projected onto the plane- (A, P4). (b) Zoom of the tridimensional region near the 
critical point. 

Fig. 7 Chaotic exit to inflation: (a) outcome of 50 orbits showed in the plane-(A, Pb) 
for Eq = 0.99999999771 and the sphere of initial conditions with radius R = 10^^. (b) 
Tridimensionl view of the region close to the critical point E for R = 10~^. (c) View of the 
plane-(A, Pb) for the case R = 10"^ and energy Eq = 0.999999999939. Note that the orbits 
remains in the neigbhorhood of E of same order of the initial sphere. 

Fig. 8 (a) Poincare surface of section Pb = 0, Pb > for A = 0.25 and Eq = 
0.92504831163435. (b) Double Poincare surface of section with the condition Pb > re- 
laxed with Eq = 0.15 and A = 2.0. It is clear the presence of destroyed regions of the phase 
space together with KAM tori structure. The tori cross the plane A = performing the so 
called cosmic cycles. 

Fig. 9 (a) Two collapsing orbits inside a torus with A = 0.25 and Eq. (b) The result 
of increasing the cosmological constant to A = 0.30, keeping the same region of the phase 
space and assuming Eq = 0.8000050001250895, is the destruction of the torus: after some 
cycles, both orbits escape to de Sitter configuration. 

Table 1 

Relation between the radius R and the gap of energy 6 Eq for the case 7 = (dust). 

Table 2 

Relation between the radius R and the gap of energy 6 Eq for the case 7 = 1/3 (radiation). 
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R 


5EI 


10-2 


0.928 X 10-4 


10-3 


0.618 X 10-6 


10-4 


0.617 X 10-8 


10-5 


0.619 X 10-1° 


10-6 


0.630 X 10-12 



Table 1 



R 


5EI 


10-2 


2.874 X 10-4 


10-3 


2.890 X 10-6 


10-4 


2.855 X 10-8 


10-5 


2.880 X 10-1° 


10-6 


2.830 X 10-12 
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